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Abstract 

A linear spin-wave approach, a variational method and exact diagonlization are used to in- 
vestigate the magnetic long-range order (LRO) of the spin-^ Heisenberg antiferromagnet on a 
two-dimensional 1 /7-depleted triangular (maple leaf) lattice consisting of triangles and hexagons 
only. This lattice has z = 5 nearest neighbors and its coordination number z is therefore between 
those of the triangular (z = 6) and the kagome {z = 4) lattices. Calculating spin-spin correlations, 
sublattice magnetization, spin stiffness, spin-wave velocity and spin gap we find that the classical 
6-sublattice LRO is strongly renormalized by quantum fluctuations, however, remains stable also 
in the quantum model. 

PACS numbers: PACS numbers 75.10.Jm, 75.50.Ee, 61.43.Hv 
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I. INTRODUCTION 



The properties of low-dimensional antiferromagnetic spin systems have been subject of 
many studies in recent years. A lot of activity in this area was stimulated by the possible 
connection of such systems with the phenomenon of high-temperature superconductivity. 
But, the rather unusual properties of quantum magnets deserve study on their own to gain a 
deeper understanding of these quantum many-body systems, especially at low temperatures. 
One of the main issues studied is the presence of long-range order (LRO) in the ground state 
of two-dimensional spin-i Heisenberg antiferromagnets (HAF), described by the Hamiltonian 

i7=j5:s,-s, (1) 

<i,j> 

on different two-dimensional lattices. The sum runs over all pairs of nearest neighbors on 
the lattice under consideration and the coupling J is positive. 

It is rather well established that LRO is present in the ground state of the spin-^ 
HAF on bipartite lattices (squared, honeycombiH, 1/5-depleted squareii, square-hexagonal- 
dodecagonali) and, contrary to some early works^'i, also on triangular latticei00. Those re- 
sults were obtained and confirmed by different methods: exact diagonalization (ED), Monte- 
Carlo simulations, spin-wave and variational approaches, series expansions and others. It 
is also worth noticing that recent experiments show that real systems can be modeled by 
spin- 1 Heisenberg antiferromagnets with different couplings on some uniformllilli and even 
depleted lattices^. 

A regular depletion of the triangular lattice by a factor of 1/4 yields the kagome lattice 
with coordination number z=4. Contrary to the triangular lattice the ground state of the 
spin-^ HAF on the kagome lattice is most likely a spin liquid. However, the kagome 
lattice is not the only regularly depleted triangular lattice. As recently has been pointed 
out by BettJll a regular depletion of the triangular lattice by a factor of 1/7 yields another 
translationally invariant lattice. The coordination number of this lattice is 2; = 5 and lies 
between those of the triangular {z = 6) and the kagome {z = 4) lattices. According to Betts 
we will call this lattice in what follows the maple leaf lattice. Since, in general, magnetic 
order is weakened by frustration and low coordination number z, it is natural to ask whether 
the magnetic LRO, present for the HAF on triangular lattice but absent for the HAF on 
kagome lattice, will survive this 1/7 depletion of the triangular lattice or not. In this paper 
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we will study this problem using several analytical and numerical methods to calculate the 
ground state of model (|I|). 

The paper is organized as follows: In Section II we briefly illustrate the geometrical 
properties of the lattice and the classical magnetic ground state, in Section III exact diago- 
nalization data for finite lattices of = 18 and 36 spins are presented and compared with 
approximate data (spin- wave and variational), in Section IV a linear spin- wave approach to 
this problem is presented, results of variational calculations are described in Section V, and 
the summary is given in Section VI. 



II. GEOMETRY OF THE LATTICE AND THE CLASSICAL GROUND STATE 

The maple leaf lattice is shown in Fig.|I[ It belongs to the class of uniform tilings in 2D 
built by a periodic array of regular polygons. In each of the equivalent sites 4 triangles and 
one hexagon meet. The maple leaf lattice has no reflection symmetry. Its unit cell (marked 
by dashed lines in Fig.^) consists of 6 sites and 15 bonds. The underlying Bravais lattice is 
a triangular one. The basis vectors are ri = h\/l a/3/2, 1/2 and V2 = [0, 1], where h is 
the distance between neigboring sites. More information can be found in Ref.0. 

The ground state of a classical spin system on such a lattice forms the starting point 
for the calculation of the ground state properties of the quantum HAF within the spin- 
wave method (section IV) and variational method (section V). As reported previously^, 
this ground state is a non-coUinear (canted) planar state with six sublattices. It can be 
characterized as follows: We denote the position of i-th hexagon (unit cell) by lattice vector 
Rj and label the sites in the unit cell by the running index n = 1, 6. Then we can write 

Si„ = s (cos (0„ + Q ■ Ri) ei sin (0„ + Q ■ Rj) eg) , (2) 

where ei and 62 are arbitrary orthogonal unit vectors. For the angles 0„ we have 
0n — 0m = ±« for nearest neighbors on the hexagon within the unit cell. The corresponding 
product of two spin vectors reads Si„, ■ S^j = cos (0„ — 0m + Q ■ (Rj ~ Rj)) • The classical 
ground state corresponds to two sets of 'wave vector' Q and pitch angle a = (f)B — 4'A, namely 
Qi = 27r^^ -j^, ^ ,ai = — fyr and Q2 = ^n^^ 0, | , 0:2 = t^r. This is a kind of trivial 
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degeneracy which one can also encounter in the system of classical spins residing on the 
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FIG. 1: 1/7-depleted triangular (maple leaf) lattice. The geometrical unit cell containing six spins 
is marked by the dashed lines. This lattice may be split into 6 equivalent triangular sublattices 
A,B,C, F. The classical ground state is represented by arrows. 

triangular lattice. The classical ground state energy is given by 

ES' = -^iV.2(l + y3), (3) 

where N is the number of sites. For spin s = \ and J = 1, we have the energy per bond 
E/hond = — (v^ + l)/20 ~ —0.137. Notice that for the HAF on triangular and kagome 
lattices E/hond = —0.125. 

The classical ground state for Q = Q2 is illustrated in Fig.|l|. One has six triangular 
sublattices A,B,...,F. The classical spins attached to the sublattices are rotated from one 
unit cell to the next one by the angle —^tt in the direction of basis vector ri and by the 
angle |7r in the direction of r2. The angle between the nearest spins on each hexagon is |7r 
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or — |7r and three spins residing on each equilateral triangle (marked by light and dark grey 
in Fig.|T]) coupled to three nearest hexagons form a 120°-structure. 

Though the classical ground state of the maple leaf lattice is more complex than that of 
the triangular lattice both are Neel states, however, with more than two sublattices. At the 
same time the classical ground state properties of the kagome lattice exhibiting a nontrivial 
ground state degeneracy are completely different. 

III. THE EXACT DIAGONALIZATION 

We used the Lanczos algorithm to calculate the lowest eigenvalue and the corresponding 
eigenstate of finite lattices of = 12, 18, 24, 30, 36 sites. This method has sucessfully applied 
to finite triangular and kagome lattices. Unfortunatly, for the maple leaf lattice only the 
lattice with A^ = 18 has the complete p6-symmetry of the infinite lattice and only multiples 
of 18 fit to the symmetry of the classical ground state. Hence we focus on the lattices with 
A^ = 18 and A^ = 36 shown in Fig.^. To reduce the Hilbert space of the Hamiltonian we 
use all possible translational and point symmetries as well as spin refiection. The number of 
symmetries of the A^ = 36 lattice is lower than that of corresponding A^ = 36 triangular and 
kagome lattices and the number of symmetrized basis in the S^^^ = ground state sector is 
378,221,361. The ground state energy per bond for A^ = 18 is Eo/bond = -0.2190042J and 
for A^ = 36 is Eo/bond = —0.2155890 J. The spin-spin correlation functions for A^ = 18 are 
collected in Table | and for A^ = 36 in Table || where they are compared to those obtained 
within spin- wave and variational approach (see below). In the fully symmetric A^ = 18 
lattice we have three different nearest-neighbor (NN) correlations. The NN correlation 
(SoSi) = —0.186299 (solid lines in Fig.^ corresponds to the classical 120° bond (see Fig.p; 
the respective averaged value for A^ = 36 is —0.177732. Both values are very close to the 
NN correlation for the triangular lattice. The NN correlation along a hexagon (dashed 
lines in Fig.^ is strong (S0S2) = —0.366673 (the respective averaged value for A^ = 36 
is —0.365555) and is close to the NN correlation of the honeycomb lattice. Finally, the 
NN correlation corresponding to a classical 90° bond (dotted lines in Fig.|^) is very small 
(SoSii) = 0.010923 (the respective averaged value for A^ = 36 is -0.037491). Hence the NN 
correlations of the quantum system reflect very well the classical ground state. 

The flnite-system order parameter corresponding to the classical ground state is the 
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structure factor (square of sublattice magnetization) 

N/a 6 



i,j=l n=l 



The values for = 18 and 36 are listed in Table A finite-size scaling of the order 
parameter with only two points seems to be not reasonable, however, doing so with a N~^^'^ 
scaling we obtain a finite value of for ^ oo. 

A better way is the direct comparison of the spin-spin correlations with those for the HAF 
on triangular and kagome latticeslll. For the presentation of the data we have to take into 
account that in the classical six-sublattice Neel state we have for instance spins with a relative 
angle of 90° leading to special correlations being zero for arbitrary distances. Therefore 
we consider as a measure for magnetic order the strongest correlations. Consequently we 
present in Fig.^ the maximal absolute correlations |(SoSr)| versus Euclidian distance r. 
As expected we have very rapidly decaying correlations for the disordered kagome case, 
whereas the correlations for the Neel ordered triangular lattice are much stronger for larger 
distances. Though the correlations for the maple leaf lattices are smaller than those of the 
triangular lattice they are significantly stronger than those of the kagome lattice and a kind 
of saturation for larger distances is suggested. 

The results for the spin-spin correlation are used to estimate the quality of the spin-wave 
and variational method used below by comparing the exact and approximate correlations 
|(SoS,)| for the finite lattice of A^ = 36 (Fig.|b). 

IV. THE LINEAR SPIN- WAVE APPROACH 

Taking into account that we have six sites in the geometrical unit cell the appropriate 
representation of the general Hamiltonian (|l]) reads 

H = J E] Sj„ ■ Sjm (5) 

<i,j\n,m> 

where i,j label the unit cells and n, m = 1, ... ,6 the different sites in one unit cell. Of 
course, the sum runs over neighboring sites, only. The linear spin-wave theory (LSWT) is 
carried out as usual. However, we need at least six different types of magnons, which makes 
the calculation more ambitious than for the triangular or the kagome lattice. We use as 
quantization axis the local orientation of the spins in the classical ground state. Performing 
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FIG. 3: The dependence of the spin-spin correlation on the Euchdean distance for the HAF 
on the finite maple leaf lattice with = 36 sites shown in Figj2|. (a) comparison between exact 
diagonahzation results for different lattices. Data for kagome are taken from Ref.i. (b) Comparison 
between exact diagonahzation, spin-wave and variational results. The lines are guide for the eyes. 
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the linear Holstein-Primakoff transformation the scalar product ■ ^jm in (§)) is replaced 
by the bosonic quadratic form 



SmSjm > S COS S COS 0^j„ ^CJ^jCtni + 

+ s (cos Q'^^ - 1 



+ s (cos e^^^ + l) (a+ + a„ia+j.) /2, 



(6) 



where n,m = 1, . . . , 6 label the different magnons in a unit cell i. represents the angle 
between the respective classical spin vectors. After transforming the Hamiltonian into 
the k-space one obtains 



H = ~Ns^ (l + v^) + Js 5] /7k, 
2 k 



(7) 



where 
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- l^k («4«ftk - a^oek) - ^73k (askOe-k - askO^k) (8) 



with 7k„ = exp(^kq„), qi = b/ V28 [-V3,5) , qa = 6/^28 (-2^3, -4^ 
6/728(373,-1) and with b being the distance between two neighboring spins. This Hamil- 
tonian can be diagonahzed by the Bogoljubov transformation 

6 

ank = ^nkOimk + ^n-k«m-k- (9) 

m=l 

The new bosonic operators a^k describe the normal modes uJmk- In order to determine them 
and the Bogoljubov coefficients one has to solve the following equations 

The solution gives six different, non-degenerated spin-wave branches — five of them are 
optical whereas the remaining one is an acoustical branch. The acoustical branch becomes 
zero in the center (k = 0) and at the edges of the Brillouin zone (k = ±Q with Q = 
^^feTf "Ts' i )■ '^^^ expansion of the zero modes in the vicinity of those points gives the 
spin-wave velocities 



^ ^7T4V39 + 23V3 
Ck=o = Jsb- 



4 (2 + x/3) 



Ck=±Q = Jsb -r . (11) 

4(7 + 4v^) 

The acoustical branch of the maple leaf lattice is similar to that of the HAF on the triangular 
lattice@'0, where one has a threefold degenerated acoustical branch being zero for k = 
and at the edges of the Brilloiun zone k = ±Q = ±27r/a l/-\/3, 1/3 . The situation for 
the HAF on the kagome lattice is completely different. Starting from the so-called classical 
k = state one obtains three branches, one dispersionsless (flat) mode tUk = and two 
degenerated acoustical branched. 

In analogy to the triangular latticed it can be shown that the zero modes k = 0, ±Q of 
the maple leaf lattice describe out-of-plane and in-plane oscillations, respectively. Therefore, 
we denote Ck=o as cy and Ck=±Q as c±. Together with the spin- wave velocity the spin 
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stiffness constitutes the fundamental parameters which determine the low-energy dynamics 
of magnetic system^. To calculate the spin stiffness p in the leading order one can 
use the hydrodynamic relation p = xc^. The magnetic susceptibilities x\\ = Xzz/V (out- 
of-plane) and x± = Xxx/V = Xyy/^ (in-plane) can be determined minimizing the classical 
energy in the limit of a vanishing external field. We find Xzz = N/ J {q + y^) , Xxx = 
Xyy = 2N/3J (^4 + y/Sj and V = 7y/3Nb'^ /12 as the volume of the lattice and from the 
hydrodynamic relation one obtains 

Pll = 0.633975Js2, = 0.211325Js2, py/p^ = 3 . (12) 

The comparison with the corresponding parameters calculated in the same order in s for 
the square and the triangular lattice are given in Table |T|. We find that the spin stiffness 
parameters for the maple leaf lattice are lower than the corresponding values of the triangular 
lattice indicating that the Neel is stronger influenced by quantum fluctuations in the maple 
leaf lattice than in the triangular one. 

The ground state energy -Eq^ is given by 

^o"^ = -^Ns (s + 1) (l + v^) + ^ ^ .;^k/2, (13) 

k m=l 

which leads in the thermodynamic limit to an energy per bond 

e^'" = (^-0.5464106s2 - 0.13652065s) J. (14) 
The sublattice magnetization 

{sDn = {S^N = ^ - If E («>nk) (15) 

k 

calculated in the thermodynamic limit is 

(S^)^ = s - 0.346. (16) 

A comparison between all those values for HAF on square, triangular and maple leaf lattices 
is given in Table |T|. Obviously, for all these parameters c, p and (S^)^ the same tendency 
is found, namely to be largest for the unfrustrated lattice and to be lowest for the frustrated 
maple leaf lattice with z = 5. Notice, that for the kagome lattice the LSWT yields diver- 
gent contributions in the sum over k in {{S'')n — s) oc J^k ('^nk'^nk) indicating a vanishing 
sublattice magnetization^!. 
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Finally, we compare the spin-spin correlations (SoSj) (where j runs over all spins in the 
system) obtained within the LSWT for the finite lattice with = 36 shown in Fig.|^ with 
the exact numerical Lanczos data (see Fig.^ and Table One finds a surprisingly good 
agreement between the approximative LSWT data and the exact Lanczos data. Hence the 
finding of finite sublattice magnetization obtained within LSWT is supported by the Lanczos 
data. 



V. THE VARIATIONAL APPROACH 



The classical ground state, described in Section II, is the basis for the construction of the 
variational Huse- ElseiS ground state which, expanded in the Ising basis states | a) of the 
total spin component = oiil, reads 

1^) = ewi^Hclass + ^Hguant + ^Hfrust) \ «)• (17) 

a ^ ^ ^ 

The operators H are diagonal in the base \a). The term 

Eclass = -%Y.'^oS] (18) 

i 

produces a proper 'classical' phase for a given state |a) in the expansion given by Eq. (p!7D. 
The sum runs over all spins in the system (i.e. the index j corresponds to a pair ii^n) of 
indices in Eqs. (^ and (|). 0j is the angle specified in Fig.|I|. 
The second operator containing variational parameters K^i^ 

Hquant = ^ Kj^SjS^ (19) 

provides the amplitude for a given basis state \a) and introduces the quantum corrections 
to the classical function by taking into account spin-spin correlations. It means that in 
this approach one starts from the state with a broken rotational symmetry and this is still 
present during the minimization procedure producing the final symmetry-broken ordered 
stateiil. 

Finally, the third operator which contains the variational parameter Ljki and the corre- 
sponding sign factors •yjki = ±1 

Hfrust = X! 'IjklLjklSjS^Sl (20) 
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describes an additional possible change of the classical phase due to the quantum fluctu- 
ations. Following the ideas of Huse and Elseiil we assume that the wave function of the 
quantum ground state (i.e. |\E') from Eq. (0) with all three terms Hdass, Hquant and Hfmst) 
has the same symmetry properties as its classical part (i.e. |\1') from Eq. (|T^ with only 
Hciass)'- the sign of the imaginary part of the wave function changes under the rotation 
-Rz(7r/3) by the angle vr/S whereas remains unchanged under the rotation Rz{2it/3) by the 
angle 2n/3 about the center of a hexagon. This transformation determines the 'shape' of 
three-spin terms Lj^i and the proper sign of •jjki = ±1 in Eq. (pO|). Similarly to the tri- 
angular latticdzl the most simple three-spin terms are 'dog legs' with j and / being nearest 
neighbors of k. For example, for spin number 5 of the = 18 lattice shown in Fig. 2 
(top) there exist four such 'interactions': -^^3,5,13, -^^3,5,4, -^^2,5,12 and -^4 5 13. Each Lij^ is con- 
nected with its corresponding 7 factor, i.e. -^3,5,13 73,5,13 = -/efa, -^^3,5,4 73,5,4 = Iefd, 
^2,5,12 72,5,12 = IBFC and 1.4,5,13 74,5,13 = 'Jdfa, where the letters A,B,C,D,E,F cor- 
respond to the 6 equivalent triangular sublattices illustrated in Fig.^. Taking into account 
that R,{n/3){ACE) = (BDF), R,{tt /3){BDF) = (CEA) and R,{2Tr/3){ACE) = (CEA), 
Rz{27t/3){BDF) = (DFB), i.e., the '120° structure' ACE (dark triangles in Fig. 1) trans- 
forms into '120° structure' BDF (grey triangles in Fig. 1) under i?^(7r/3) or into itself under 
Rz{27t/3), one obtains a proper signs of 7 factors 

laa/S = -Ippa- (21) 

The index aa(3 means that two (different) spins in the three-spin term belong to the same 
'120° structure', the remaining one belongs to the other '120° structure'. Thus, for example, 
if one puts 7_bfc = 1 it follows that 7cad = —1 (or 72,5,12 = 1 and 712,13,10 = —I, see Fig. 
2). 

How does one choose the variational parameters Kjk and Ljki in Eq. (0) for the HAF 
on the maple leaf lattice? We have applied two criteria: a better choice of parameter 
space should give a lower value of the ground state energy and, if two energies for different 
parameter spaces are approximately the same, one should choose the parameter space which 
leads a lower value of the variance {H'^) — (H)'^. In order to find an optimal choice of 
the wave function we have tested some possibilities for the parameter space for the fully 
symmetrical A^ = 18 lattice taking into account the whole = basis in the expansion 
(P!7|). The best choice found is the following five-parameter space (results for the correlation 
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FIG. 4: Variational energy per bond for the spin system on the 1/7-depleted triangular (maple 
leaf) lattice as a function of N~^^'^. 

are collected in Table I): {Khex, Ktr, Kothers, Li, L2). Spins 'interacting' via K^ex are nearest 
neighbors lying on a hexagon: hex = AB, BC, CD, FA, those 'interacting' via Ktr are 
nearest neighbors belonging to the 120° structure: tr = EC, CA, AE or BF, FD, DB; and 
all remaining nearest neighbors are coupled by Kothers, thus others = BE, BC, EF. Note 
that there is no long-range variational parameter for pairs of spins not being the nearest 
neighbors. Moreover, one takes into account only three, from the four existing 'dog leg' 
interactions, i.e., 'dog legs' around a hexagon are absent. For example, in each point F one 
has Li = Lbfc and L2 = Lefd = Lafd and Lefa is absent (or correpondingly for spin 
number 5 in Fig. 2, Li = ^2,5,12 and L2 = -^^3,5,4 = 1^13,5,4 and ^3,5, 13 is absent). All the 
expectation values of operators reported in the following are calculated for this choice of the 
variational parameters. 

Having obtained the ground state function one can calculate the expectation values of 
the operators which characterize the ground state of a given, finite spin system. This can 
be accomplished by a Monte- Carlo approach^ and the finite size scalingil tells how to 
extrapolate those expectation values to the thermodynamic limit. We have investigated the 
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FIG. 5: Square of sublattice magnetization, function of N Circles - values obtained 

by applying the variational method. Solid line - fit to the circles. Sizes of circles are comparable 
to the statistical error bars. 

finite systems of 18, 72, 162, 288 spins with periodic boundary conditions. Note that they 
have the full symmetry of the maple leaf lattice. The relevant quantities are collected in 
Table |I^ and the finite size analysis is presented in Figs.§ and ^. 

The leading term of the finite-size correction of the ground state energy per bond e is 
The data in Table can be fitted to this dependence (see Fig.^ and hence the 
energy per bond Coo in the thermodynamic limit is obtained: e(A^) = Coo + aN^'^^"^ with 
Coo = — 0.1988(2)J and a = —0.7327(164) J. This value for e^o is about 3% higher than the 
value obtained from spin- wave theory (see Eq. (0)). 

In Fig.^ the finite-size extrapolation of the square of sublattice magnetization defined 
in Eq. (g) is shown. We find m^{N) = + cA^-^/^ + dN-^ with = 0.0723(10), 
c = 0.0444(18) and d = 3.650(60) suggesting that the long-range magnetic order persists 
in the ground state of this spin system. Note, however, that the applied variational ansatz 
tends to overestimate the magnetic order (see Ref.il and table |l| as well as Fig.^a). 

The variational approach enables us to calculate the spin gap A = Eq — Ei, where Eq 
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FIG. 6: The spin-gap A = Ei — Eq vs. in units of J. Eq and Ei denote variational energies 
of the ground state and a first excited state, respectively. Errors result from adding the errors for 
Eq and El. 

{El) is the variational energy in the subspace of total Sz = {Sz = 1). This new aspect 
of the Huse-Elser ansatz was used for the first time in Ref.i to calculate the spin gap for 
the HAF on the square-hexagonal-dodecagonal lattice. Magnetic LRO is connected with 
gapless Goldstone modes whereas quantum disorder in the ground state is accompanied 
by a finite spin gap. Therefore the calculation of A yields an additional argument for or 
against the existence of magnetic LRO order in the ground state. Fig.|^ shows the finite- 
size extrapolation of the spin gap according to the relation A(A^) = Aqo + dN"^ with 
Aoo = —0.0180(10) J and d = 6.3610(34) J. The negative Aqo is a result of the limitted 
accuracy of the approximation but nevertheless suggests a zero spin gap. Hence we have an 
additional indication for the existence of LRO. 
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VI. SUMMARY 



In this paper the results of exact diagonahzation, hnear spin-wave theory and a Huse-Elser 
hke variational investigation for the ground state of the spin-| Heisenberg antiferromagnet on 
a new 1/7-depleted triangular (maple leaf) lattice are presented. The coordination number 
of this frustrated lattice z = 5 lies between those of the triangular and the kagome lattices. 
Quantum fluctuations and frustration tend to destroy classical magnetic ordering. Their 
influence becomes the stronger the smaller the coordination number. But contrary to the 
kagome lattice with z — 4 ior the maple leaf lattice we find strong arguments that the 
classical six-sublattice Neel LRO survives the strong quantum fiuctuations present in this 
frustrated quantum magnet. This conclusion is drawn from the calculated values of the 
spin-spin correlation, sublattice magnetization, spin stiffness, spin-wave velocity as well as 
the spin gap. 

The comparison between exact data and approximate data for the spin-spin correlation on 
finite lattices gives a surprising well agreement between the linear spin-wave and the exact- 
digonalization data whereas the variational approach tends to overestimate the strength of 
correlations. 

Finally, we mention that on the passage from the triangular to the 1/7 depleted (maple 
leaf) lattice (i.e., some interactions J in spin system on triangular lattice are varied 
from J = 1 to J = 0), one would encounter a transition between thrce-sublattice and 
six-sublattice Neel LRO which may have interesting features worth to be considered in future. 
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TABLE I: The exact values of the spin-spin correlation (SoSj) compared to their variational values 
for the lattice of = 18 sites (see Fig.Q). 



j exact variational j exact variational 

0.750000 0.750000 7 0.180027 0.200775 

1 -0.186299 -0.180068 8 0.140873 0.162808 

2 -0.366673 -0.343444 9 -0.072868 -0.106613 

4 0.039003 0.021877 11 0.010923 0.005425 

5 0.145098 0.171218 17 -0.174804 -0.183760 

6 -0.099672 -0.106613 
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TABLE II: The exact values of the spin-spin correlation (SoSj) compared to their spin-wave and 
variational values for the lattice of = 36 sites (see Fig.|2|). Statistical errors are given in paren- 
theses. 



3 


exact 


spin-wave 


variational 


j 


exact 


spin-wave 


variational 


1 


-0.1154 


-0.16759 


-0.1681(90) 


19 


0.0660 


0.05816 


0.1574(70) 


2 


-0.3418 


-0.31894 


-0.3408(90) 


20 


0.1458 


0.14131 


0.1603(60) 


3 


-0.2008 


-0.19143 


-0.1703(90) 


21 


-0.0111 


-0.00801 


-0.0784(90) 


4 


0.0618 


0.03135 


0.0249(70) 


22 


-0.3929 


-0.31878 


-0.3395(90) 


5 


0.1394 


0.15200 


0.1701(50) 


23 


-0.0433 


-0.06302 


0.0027(60) 


6 


-0.0155 


-0.00625 


-0.0811(90) 


24 


0.0434 


0.01186 


0.1503(70) 


7 


-0.0243 


-0.01525 


-0.0788(90) 


25 


-0.0491 


-0.03622 


-0.0845(90) 


8 


0.0089 


-0.00304 


0.0190(60) 


26 


-0.0448 


-0.02126 


-0.1393(90) 


10 


0.1142 


0.12745 


0.1546(60) 


28 


0.0034 


-0.01380 


0.0089(60) 


11 


-0.0493 


-0.04012 


-0.1470(90) 


29 


0.0298 


0.01152 


0.1249(70) 


12 


-0.0155 


-0.00625 


-0.0832(70) 


30 


-0.1059 


-0.09804 


-0.0990(90) 


13 


0.1488 


0.12892 


0.1960(50) 


31 


-0.0740 


-0.06787 


-0.0939(80) 


14 


0.0327 


0.01839 


0.1261(80) 


32 


0.0387 


0.02060 


0.0156(70) 


15 


-0.0797 


-0.06993 


-0.0946(80) 


33 


0.1785 


0.15290 


0.1997(50) 


16 


-0.0500 


-0.03038 


-0.1407(90) 


34 


0.0501 


0.04121 


0.1313(70) 


17 


0.0390 


0.02791 


0.0151(60) 


35 


-0.1561 


-0.11596 


-0.1716(90) 


18 


-0.1059 


-0.09804 


-0.0992(90) 










TABLE III: Comparison of the LSWT results for the spin-wave velocities, spin stiffness parameters 


and sublattice magnetization for 


the squareHH, 


the triangularis 


E and the maple 


leaf lattice 


(J = l 


,6=1 and 


s = l/2). 












lattice 


C|| 


c± 




P\\ 


Pi. 




square 


1.4142135 


1.4142135 




0.25 


0.25 


0.304 


triang 


ular 


1.2990381 


0.9185586 




0.2165063 


0.1082532 


0.239 


maple 


leaf 


1.1127356 


0.6774616 




0.1584936 


0.0528312 


0.154 
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TABLE IV: The ground state energy per bond Eo/bond, the square of sublattice magnetization 
and the spin gap for the HAF on finite 1/7-depleted triangular (maple leaf) lattices (J = 1). For 
the N = 18 and the iV = 36 lattice the results of exact diagonalization are also included. In the 
case of the N = 18 lattice the variational values were obtained in the whole basis of Ising states, 
for larger systems the Monte-Carlo method was applied. Statistical errors are given in parentheses. 



N 




Eg/bond 


2 

m 


gap 


18 


exact 


-0.2190 


0.2832 


0.5452 




variational 


-0.2083 


0.2855 


0.3353 


36 


exact 


-0.2155 


0.1534 






variational 


-0.2027(1) 


0.179(1) 


0.146(5) 


72 




-0.2001(1) 


0.128(1) 


0.071(7) 


162 




-0.1991(1) 


0.099(1) 


0.020(10) 


288 




-0.1988(1) 


0.088(1) 


0.005(12) 


GO 




-0.1988(2) 


0.072(1) 


-0.019(25) 
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